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The  work  reported  in  this  note  was  done  for  the  purpose  of  determining 
the  feasibility  of  neutralizing  the  potential  difference  that  occurs  on  a 
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by  emitting  ions  from  the  positively  charged  surface  region  of  the  spacecraft 
and  by  letting  this  current  impact  on  the  negatively  charged  surface  region. 
In  order  to  Investigate  the  main  features  of  the  proposed  scheme,  we  made 
some  simplifying  assumptions  about  the  geometry  of  the  problem.  We 
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I.  INTRODUCTION 


The  work  reported  in  this  note  was  done  for  the  purpose  of  determining  the 
feasibility  of  neutralizing  the  potential  difference  that  occurs  on  a  spacecraft. 

It  has  been  proposed  that  neutralization  can  be  accomplished  by  emitting  ions 
from  the  positively  charged  surface  region  of  the  spacecraft  and  by  letting  this 
current  impact  on  the  negatively  charged  surface  region.  In  order  to  investigate 
the  main  features  of  the  proposed  scheme,  we  made  some  simplifying  assumptions  about 
the  geometry  of  the  problem.  We  approximated  the  geometry  by  two  conducting 
infinite  half  planes  and  investigated  the  trajectories  of  a  positive  ion  emmitted 
at  some  point  on  the  positively  charged  infinite  half  plane  with  given  initial 
velocity  and  a  given  direction  of  emission. 

II.  MATHEMATICAL  STATEMENT  OF  THE  PROBLEM 


We  consider  the  plane  y=o  and  assume  a  cut  on  this  plane  along  the  z  axis. 
Let  there  be  a  given  potential  difference  V  between  the  two  half  planes  x  >  o 
and  x  <•  o. 

The  equations  of  motion  of  a  charged  particle  in  the  electric  field  produced 
in  such  a  configuration  are 


mdF  = 


(1) 
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where  the  initial  velocities  of  the  particle  are  v  and  v  ,  and  the  initial 
position  is  x  =  xQ,  y  =  o. 

In  polar  coordinates 
x  =  r  cos  G 
y  *  r  sin  0 

the  equations  of  motion  become 

ft 

(¥0  -  rQ2)  cos  e  -  (2re  +  r8)  sin  0  -  -  —  sin  0 

(2r0  +  r0  )  cos  0  +  (r  -  r02 )  sin  0  =  +  ^  cos  ®  (3) 


with 


mir 

These  equations  can  be  combined  into 

*  *  "  ft 

2  r0  +  r0  =  -- 

r 

and 

r  -  r62  =  0 

Eq.  (5)  can  be  simplified  to 


(4) 

(5) 
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|t(r2e  )  =  8  (7) 

which  provides  the  angular  momemtum  integral 

r2  0  =  r  v  +  8  1  (8) 

o  yo 

Eq.  (6)  shows  that  r  >  0  always,  thus  if  r  (0)  _>  0,  r  itself  also  increases  all  the 
time. 

The  energy  integral  is  obtained  from (6)  through  multiplication  by  r,  and  using  (7): 

r2  +  r2  62  )  -  6  6  =  k  v  2  (9) 

o 

V2  =  V2  +  V2  , 

o  yo 

r  =  x 

o  o 

Eliminating  9 

equation  for  r: 

d2r 
dtr 

which  becomes 

d2f 

d^7 

after  making  the  transformation 

u  *  8t  +  r  v 

o  yo 

f  -  /IT  r  (12) 


from  (6)  with  the  aid  of  (8),  we  obtain  a  second  order  differential 


(Bt  +  r0vyo ) ‘ 
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A  second  order  differential  equation  for  the  angle  9 
although  it  is  of  limited  interest  only: 


.d0  d29.2 

^du  "  U  du2^ 
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can  also  be  obtained, 


(13) 
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The  differential  equation  of  the  trajectories  can  be  obtained  by  eliminating 
t  (or  u) .  From  (7)  we  have 


r29 
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(14) 


Eliminating  8  with  the  aid  of  (8),  this  results  in 

2 

UV„ 
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III.  NUMERICAL  SOLUTIONS 


The  equations  of  motion  (1)  can  be  cast  in  a  dimensionless  form 
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y 

x2  +  y2 
x 

x2  +  y2 


with  the  initial  conditions  that  at  time  t=o 


x  =  1 

y  =  0 


V 

X 

=3 
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xo 

V 
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f  »  ev 

where  the  initial  velocities  v  and  v  are  in  units  of  v  8  *  6  =  —  • 

XO  yo  IWT 

We  solved  numerically  this  set  of  differential  equations  with  the  help  of  a 
subroutine  DHPCG  based  on  the  Hamming's  Modified  Predictor  Corrector  Method.  We 

assumed  various  initial  conditions  for  v  ,  v  . 

xo  yo 

In  the  Table  we  give  some  results  for  cases  where  there  is  no  initial  velocity 
in  the  perpendicular  direction,  that  is  v  -  0,  and  the  inital  horizontal 


a r 
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motion  is  towards  the  center.  Also  v  E  0.  x.  is  the  distance,  from 

zo  * 

the  center,  of  the  point  where  the  positive  ion  impacts  on  the  negatively 

charged  half  plane,  y  is  the  maximum  value  of  the  y  coordinate  of  the 

max 

particle  along  its  trajectory.  Note  that  our  results,  up  to  -v  'v  3.5  are 

xo 

smooth  in  the  sense  that  as  -v  increases,  x  .  decreases.  At  -v  =  4, 

xo  *  xo 

and  higher  values,  the  results  are  not  reliable.  Whether  the  problem  is  in  the 

numerical  procedure  is  not  clear  at  the  present  time. 

Fig.-s  1,  2,  and  3  portray  families  of  particle  trajectories  with  different 

initial  conditions.  In  Fig.  1,  the  initial  velocity  is  in  the  x-direction.  It 

is  clear  that  as  v^q  increases,  the  impact  distance, x^  is  drastically  reduced. 

Whether  there  is  an  optimum  initial  velocity,  corresponding  to  a  minimum 

impact  distance,  cannot  be  determined  before  the  computational  problem  referred 

to  above,  is  clarified.  In  Fig.  2,  the  initial  velocity  is  in  the  v-direction. 

Here  the  impact  distance  increases  with  increasing  velocity.  In  Fig.  3,  the 

absolute  value  of  the  initial  velocity  is  kept  constant,  while  the  direction  of 

the  emission  is  varied.  There  is  a  critical  angle,  0  ,  at  which  the  impact 

o 

distance  is  the  same  as  for  v  =0.  For  9  <  0  the  impact  distance  is  reduced, 

o  o 

for  0  >  9  it  is  increased.  With  the  chosen  v  =  0.1,  0  'v  83°.  Fics.  4,  5. 

o  o  o—  "  ’  * 

6  give  a  more  detailed  picture  of  the  dependence  of  the  impact  distance  and  of 
the  maximum  height  of  the  trajectory  on  the  parameters. 
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CAPTIONS  FOR  FIGURES 

FIG.  1  Trajectories  for  various  initial  velocities  in  the  horizontal 

direction.  The  impact  distance  decreases  with  increasing  velocity 

FIG.  2  Trajectories  for  various  initial  velocities  in  the  vertical 

direction.  The  impact  distance  increases  with  increasing  velocity 


FIG.  3  Trajectories  for  various  initial  emission  angles.  The  impact 
distance  increases  with  increasing  angle. 

FIG.  4  The  impact  distance  (x^)  and  the  maximum  height  of  the  trajectarv 
(ym)  as  functions  of  the  initial  vertical  velocity. 

FIG.  5  The  impact  distance  (x^)  and  the  maximum  height  of  the  trajectary 
(ym)  as  functions  of  the  initial  horizontal  velocity. 

FIG.  6  The  impact  distance  (x^)  and  the  maximum  height  of  the  trajectary 
(ym)  as  functions  of  the  angle  of  emission.  Note  that  there  is 
a  critical  angle  0  ,  where  x.  (60,  v0)  =  x^  (vQ  =  0). 


FIG.  5 


